Quantum speed limit for non-Markovian dynamics 
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We derive a Margolus-Levitin type bound on the minimal evolution time of an arbitrarily driven 
open quantum system. We express this quantum speed limit time in terms of the operator norm of 
the nonunitary generator of the dynamics. We apply these results to the damped Jaynes-Cummings 
model and demonstrate that the corresponding bound is tight. We further show that non-Markovian 
effects can speed up quantum evolution and therefore lead to a smaller quantum speed limit time. 
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What is the maximal speed of evolution of a quan- 
tum system? This question is of fundamental importance 
in virtually all areas of quantum physics, ranging from 
quantum communication pQ , computation 2J and metrol- 
ogy [3] to optimal control theory [3] and nonequilibrium 
thermodynamics [5]. For closed systems, the quantum 
speed limit time which determines the maximum rate 
of evolution can be obtained by combining the results of 
Mandelstam-Tamm (MT) [5] and Margolus-Levitin (ML) 
0: it is given by t QS l = ma,x{nh/{2AE), nh/(2E)}, 
where AE is the variance of the energy of the initial state 
and E its mean energy with respect to the ground state. 
The minimal time a quantum system needs to evolve from 
a given state to an orthogonal state is thus determined by 
its initial energy [8HT5]. The latter reflects the fact that 
the Hamiltonian is the generator of unitary Schrodinger 
dynamics. It is worth emphasizing that the existence of 
a speed limit time is a purely quantum effect which van- 
ishes when h goes to zero. Generalizations of the MT 
and ML findings to nonorthogonal states and to driven 
systems have been provided in Refs. [TlHTrj] and [P7HT§] . 
Recently, the quantum speed limit time has been derived 
for open systems described by positive nonunitary maps; 
applications to dephasing in noisy channels and quan- 
tum parameter estimation have been discussed [5D1 |5T] . 
In both approaches, the speed limit time was obtained in 
terms of the variance of the generator of the evolution, 
which reduces to the MT expression in the case of closed 
system dynamics. To our knowledge, no ML type bound 
has been proposed for open quantum systems to date. 

In this paper, we use a geometric approach to derive a 
quantum speed limit time valid for open system dynam- 
ics with possibly time-dependent nonunitary generators. 
In contrast to previous studies, we obtain a bound that 
depends on the mean of the generator; it therefore re- 
duces to the ML formula for unitary processes. MT type 
bounds are usually derived with the help of the Cauchy- 
Schwarz inequality. The latter invariably leads to ex- 
pressions containing the variance of the generator, and 
is hence not suitable for getting ML type bounds. We 
here solve this technical challenge by making use of the 



von Neumann trace equality [22- 24j . In the following, 
we obtain a quantum speed limit time that depends on 
the operator norm of the nonunitary generator and show 
that the so obtained ML bound is not only sharper than 
the MT bound, it is also tight. By employing both in- 
equalities, we are able to obtain a unified quantum speed 
limit time for generic positive open system dynamics. We 
further apply these results to investigate the influence 
of non-Markovianity on the rate of quantum evolution. 
Non-Markovian (or memory) effects become important 
when the relaxation time of the system is comparable to 
the relaxation time of the environment [25l [26]. They 
have been shown to play a central role in the creation 
of steady state entanglement [27J and in the description 
of quantum coherence in photosynthetic systems [25]. A 
recent experiment with photons in a controllable non- 
Markovian environment has been reported in Ref. [29j . 
Interestingly, we will show that non-Markovian dynam- 
ics can lead to smaller quantum speed limit times. 

Geometric approach. We consider a possibly driven 
open quantum system and ask for the minimal time that 
is necessary for it to evolve from an initial state pa to a 
final state p T . Without loss of generality, we assume that 
the initial state is pure, po — |?/'o)(V'o|- The case of an 
initially mixed state can be treated by purification in a 
sufficiently enlarged Hilbert space [3U]. Note that under 
nonunitary dynamics, the final state p T will be generally 
mixed. The basis of our geometric approach is provided 
by the Bures angle £(po,p T ) between initial and final 
states of the quantum system [301 GH] > 

£(Po,Pt) = arccos \J (V'oIPtIV'o)) • (1) 

The Bures angle is a generalization to mixed states of the 
angle in Hilbert space between two state vectors [52"] . 

To evaluate the quantum speed limit time, we consider 
the dynamical velocity with which the density operator 
of the system evolves [TUJ . The latter is given by the time 
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derivative of the geometric Bures angle ([TJ) , 



df 



C(fto,pt) < 
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(2) 



Using the definition (fTl), Eq. ([2]) can be written as 

2cos(£)sin(£)£< KVolPtlV'o}!- (3) 

Expression ([3]) will serve as the starting point for our uni- 
fied derivation of ML and MT type bounds on the rate of 
quantum evolution, using respectively the von Neumann 
trace inequality and the Cauchy-Schwarz inequality. 

Margolus-Levitin bound. To illustrate the use of the 
von Neumann trace inequality, we begin by providing a 
derivation of the ML bound in the case of driven uni- 
tary dynamics. The generator is here given by the time- 
dependent Hamiltonian H t of the system, and its density 
operator obeys the von Neumann equation, 
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Substituting Eq. Q into Eq. we have, 

2 cos (£) sin {C)t<\\ | [H t , p t ] \ip )\- 

Equation ^ can be estimated from above with the help 
of the triangle inequality and we obtain, 



(5) 



2cos(£)sin(£)£ < - (|tr {H t p t p }\ 



\ix{ptH t po}\). 



(6) 

To proceed, we introduce the von Neumann trace in- 
equality for operators which reads 



(7) 



|tr{AiA 2 }| < y^at^i 



=i 



Inequality ([T]) holds for any complex n x n matrices A\ 
and Ai with descending singular values, o~\.\ > ... > <7\. n 
and (T2.1 > ••• > <?2,n- The singular values of an operator 
A are defined as the eigenvalues of vA^A [33]. For a 
Hermitian operator, they are given by the absolute value 
of the eigenvalues of A, and are positive real numbers. 
If Ai and A-i are simple (positive) functions of density 
operators acting on the same Hilbert space, Eq. Q re- 
mains true for arbitrary dimensions [24] . The singular 
values of the operators A and A> are moreover identical 
By taking A = {H t p t ) ] 
Eqs. ^ and Q, we thus find, 



p\h\ = ptH t , and combining 



2 cos (£) sin (£) £ < — } 
h t—* 1 



o-i, 



(8) 



where <7j are the singular values of H t p t and Pi = 5n 
those of the initial pure state po- For a Hermitian oper- 
ator, the operator norm is given by the largest singular 



value, milop = 01, while the trace norm is equal to their 
sum, ||A|| tr = ^Zi^i 33 . As a consequence, 



2cos(£)sin(£)£ < - \\{H t p t ) 
n 



|o P <^||(ff t pt)|| tr . (9) 



We note, in addition, that the trace norm is given by, 

\\(H t p t )\\ tI = ti{\H t p t \} = (H t ), (10) 

when the instantaneous eigenvalues of H t are all positive 
(the latter can always be realized by properly choosing 
the zero of energy [7]). Integrating Eq. ([£]) over time from 
t = to t = t, we arrive at the inequality, 



T > 



2E T 



sin 2 (£(p,p T )), 



(11) 



with the time-averaged energy E T — (1/r) J Q dt (H t ). 
Equation (111 is the ML bound for driven closed systems. 



The above derivation can be easily extended to arbi- 
trary time-dependent nonunitary equations of the form, 



Pt = L t {p t ) . 



(12) 



with positive generator L t [33] . The latter are trace class 
(super-)operators in a complex Banach space and need 
generally not be Hermitian. However, for symmetric 
norms, such as the Schatten p-norm, ||£t|| P = Ei ^i] 1 ^: 
that we here consider [31], ||£J|| — 1 1 -Ml- As a result, all 
previously used definitions and inequalities remain valid 
[35] , Substituting Eq. ([TJ) into Eq. ^jj, we then have, 



2cos(£)sin(£)£< \(^ \L t (p t )\^o) 



(13) 



which is the nonunitary generalization of Eq. ([5| . Noting 
that (if) \L t (p t )\if) ) = tr{L t (/9 t )p } and employing the 
von Neumann trace inequality Q, we obtain, 



2 cos (£) sin (£) £ < \_. ^iPi = ^i i 



(14) 



where Ai are the singular values of the operator L t (pt 



Equation ( 14 ) can again be estimated from above by the 

,. (15) 



operator norm and the trace norm to yield, 



2cos(£)sin(£)£ < \\(L t ( Pt 



< 



(LtiPt 



Integrating Eq. (|15|) over time, we eventually find, 
t > max<! ~t~hk , -T^rr )■ sin 2 (£(p, p T )), 



A? ' A* 



(16) 



where we have defined A° p ' tr = (1/t) dt ||£t(pt)|| p,tr- 
Equation ( 16 ) provides a ML type bound on the rate 



of quantum evolution valid for arbitrary positive driven 
open system dynamics. 

Mandelstam-Tamm bound. We next derive a unified 
bound for the quantum speed limit time for open systems 



3 



by generalizing the method presented in Ref. [2L based 
on the relative purity. To this end, we rewrite Eq. ([3| as 



2cos(£)sin(£)£ < |tr {L t (p t ) p }\ 



(17) 



The latter can be estimated from above with the help of 
the Cauchy-Schwarz inequality for operators: 



2 cos (£) sin (£) £ < yj tr {L t (p t ) L t {p t )^} tr{pl). (18) 
Since po is a pure state, tr {p 2 ,} = 1, an d we obtain, 



2cos(£)sin(£)£ < ^tv{L t {p t ) L t (p t )t} = ||it(pt)|| hs , 

(19) 

where ||A|| h8 = ^trjAtA} = gf is the Hilbert- 

Schmidt norm |33| . Integrating Eq. ( |11[ ) over time leads 
to the following MT type bound for nonunitary dynamics, 



t > ^sin 2 (£(p,p T )), 



(20) 



where A^ s = (1/r) di ||L t (p t )||h s . For unitary pro- 
cesses, A^ s is equal to the time-averaged variance of the 
energy. For initially pure states, relative purity and fi- 



delity are identical, and the bound (20) thus reduces to 
the one derived in Rcf. 21 , since sura; > |cosa; — 1|, 
< x < 7r/2; they are, however, different for initially 
mixed states. Combining Eqs. (16) and (20), we obtain, 



T QSL 
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1 1 

A^' A*? 



1 

aT s 



sin 2 (£(p,p r )). (21) 



Equation (21 ) provides a unified expression for the quan- 



tum speed limit time for generic (positive) open system 
dynamics. It represents a general extension of the MT 
and ML result, r QSL = ma,x{nh/(2AE),Trh/(2E)}. 

We may go a step further by noting that for trace class 
operators the following inequality holds (see Ref. [35] . 
Theorem 1.16), 



< 



lis 



< 



(22) 



As a result, 1/A° P > 1/A^ S > 1/A* r , and we can therefore 
conclude that the ML type bound based on the operator 
norm of the nonunitary generator provides the sharpest 
bound on the quantum speed limit time. We will show 
below that the bound can be attained and is hence tight. 

Non-Markovian effects. We may use the above results 
to investigate the influence of non-Markovian dynam- 
ics on the quantum speed limit time. To this end, we 
consider the exactly solvable damped Jaynes-Cummings 
model for a two-level system resonantly coupled to a 
leaky single mode cavity |36[ I37j ; the environment is sup- 
posed to be initially in a vacuum state. The nonunitary 
generator of the reduced dynamics of the system is 
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1 



Lt(pt) = It cr^p t CT + - - o+o- p t - -p t cr+cr 



(23) 



where a± — a x ± ia y are the Pauli operators and 7t 
the time-dependent decay rate. By assuming that there 
is only one excitation in the combined atom-cavity sys- 
tem, the environment can be described by an effective 
Lorentzian spectral density of the form, 



i 



7o A 



2tt (w - w) 2 + A 2 ' 



(24) 



where ujq denotes the frequency of the two-level system, 
A the spectral width and 70 the coupling strength. The 
time-dependent decay rate is then explicitly given by, 



It 



27 Asinh(dt/2) 
dcosh(eft/2) + Asinh(dt/2) ' 



(25) 



where d = -\/A 2 — 270A. In the interaction picture, the 
reduced density operator of the system at time t reads, 



Pii(0) 



Pt = 



- /o dt 'lt> 



p 10 (0)e-/o rft 'W2 



Pt (0)e-^V/ 2 l-pn(0) 



■ Jo dt' lt , 



(26) 

We shall examine the case where the system starts in the 
excited state, pn(0) = 1 and pio(0) = 0. For vanishing 
coupling, the system is isolated and in a stationary state. 
For finite coupling, the two-level system is driven by the 
bath. The correlation time of the bath is tb = A -1 , while 
the decay time of the system is equal to = 7Q -1 - The 
non-Markovian properties of the model have been inves- 
tigated in Refs. [55H4T] : The dynamics is Markovian in 
the weak-coupling regime, 70 < A/2. For large time scale 
separation, tb <C T5, or equivalcntly 70 <C A, the de- 
cay rate is constant, j t — 79. The dynamics becomes 
non-Markovian for strong coupling, 70 > A/2, which cor- 
responds to an imaginary parameter d. In this regime, 
the decay rate is an oscillatory function of time. 



Figure 1(a) shows the quantum speed limit time (21) 
for the two-level system, as a function of the coupling 
strength 70, obtained for the three different norms, in 
the case r = 1. We can distinguish two different phases. 
The speed limit time exhibits a plateau independent of 70 
for moderate coupling and then decreases for large cou- 
pling amplitudes. Our second observation is that the ML 
bound based on the operator norm is sharper than the 
MT bound based on the Hilbert-Schmidt norm, in agree- 



ment with Eq. (22). It is also sharper that the ML bound 



based on the trace norm. Remarkably, the operator-norm 
bound is tight as it reaches the actual driving time r over 
a large range of coupling strengths. 

The above behavior can be explained by evaluating the 
singular values of the operator p T in the strong Marko- 
vian limit 7 T = 70 . The two values are equal and given by 
I Pill = |poo| = l7rexp(- J r dtj t )\ = 7oexp(-7o"r). For 
small coupling, 70T <ti 1, the singular values are thus pro- 
portional to the coupling strength, \pn \ ~ 70- For larger 
coupling such that 70T 3> 1, the singular values are in- 
dependent of 70, |pn| — 0. The plateau seen in Fig. 1 
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FIG. 1: (color online) Quantum speed limit time tqsl, 



Eq. (21 1, for the damped Jaynes-Cummings model as a func- 
tion of the coupling strength. The three curves are obtained 
for the trace norm (green dotted), Hilbert-Schmidt norm (blue 
dashed) and the operator norm (red solid) of the nonunitary 
generator L(p t ), Eq. ( |23[ ). Parameters are A = 50, luq — 1 
and t — 1. 
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FIG. 2: (color online) Time-averaged operator norm of the 
nonuni tary generator L(p t ) of the damped Jaynes-Cummings 
model (23 1. The dotted line corresponds to the Markovian 
The inset shows the fidelity F(po,p T ) = 
Same parameters as in Fig. 1. 



result, Eq. 
cos£(po, p 
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is hence a signature of Markovian dynamics and follows 
from the time independence of the decay rate. The height 
of the plateau can be determined by computing the time 
averaged norm of p\ plotted in Fig. 2, 



dt\\p t \ 



[1 - exp(-7 r)] , 



(27) 



with n = 1, V2 and 2 for the operator, Hilbert-Schmidt 
and trace norms, respectively. The constant n is equal 
to 2 1 /?' for the general Schatten p-norm. Equation (|27|) 
hence shows that the operator norm (p = oo) is the only 
p-norm for which the plateau reaches the actual driving 
time r. Furthermore, the increase of the norm of the rate 
p t in the strong coupling regime, 70 > A/2, appears as a 



consequence of the (oscillatory) time dependence of the 
decay rate 74, and is thus a purely non-Markovian effect. 
We therefore reach the interesting conclusion that non- 
Markovian dynamics can increase the rate of evolution of 
a quantum system, and thus reduce the quantum speed 
limit time below its Markovian value. 

Conclusions. We have derived a quantum speed limit 
time that generalizes the familiar MT and ML results 
to generic time-dependent (positive) dynamics of open 
quantum systems. In particular, using the von Neu- 
mann trace inequality, we have obtained an expression 
of the speed limit time in terms of the operator norm 
of the nonunitary generator of the evolution. We have 
demonstrated that the latter bound is sharper than any 
bound based on a Schatten p-norm, such as the trace and 
Hilbert-Schmidt norms, and that it is moreover tight. 
Applying these results to the damped Jaynes-Cummings 
model has additionally shown that non-Markovian ef- 
fects can lead to faster quantum evolution, and hence 
to smaller quantum speed limit times. 
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